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This chapter further extends the results obtained in chapters 4 and 5 (from linear 
equation to linear systems). Each algorithm is thoroughly proved and then an example is 
given. 

Five integer number algorithms to solve linear systems are further given. 



Algorithm 1 . (Method of substitution) 

(Although simple, this algorithm requires complex computations but is, 
nevertheless, easy to implement into a computer program). 

Some integer number equation are initially solved (which is usually simpler) by 
means of one of the algorithms 4 or 5. (If there is an equation of the system which does 
not have integer systems, then the integer system does not have integer systems, then 
Stop.) The general integer solution of the equation will depend on n - 1 integer number 
parameters (see [5]): 



(A) 




i = 1,// , 



where all functions f- U are linear and have integer number coefficients. 



This general integer number system (p { ) is introduced into the other m — 1 
equations of the system. We obtain a new system of m - 1 equations with n - 1 unknown 
variables: 



k{ l) , h = %n-\, 

which is also to be solved with integer numbers. The procedure is similar. Solving a new 
equation, we obtain its general integer solution: 

(A) = 4 2) ( /t l (2) ’-’^-2)’ * 2 = 

where all functions /)'; 2 ' are linear, with integer number coefficients. (If, along this 



algorithm we come across an equation, which does not have integer solutions, then the 
initial system does not have integer solution. Stop.) 

In the case that all solved equations had integer system at step (j), 1 < j <r 
( r being of the same rank as the matrix associated to the system) then: 



(Pj) -1) = (k\ [ j) ’~’ k n-j ), ij = U —j + 1 , 



/ i ' < / 1 are linear functions and have integer number coefficients. 
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Finally, after r steps, and if all solved equations had integer solutions, we obtain 
the integer solution of one equation with n — r + l unknown variables. 

The system will have integer solutions if and only if in this last equation has 
integer solutions. 

If it does, let its general integer solution be: 

(P r ) k^~ l) =4 <r) k\ r \...XZ ,i,. = l,n-r + l, 
where all f. r) are linear functions with integer number coefficients. 

We’ll present now the reverse procedure as follows. 

We introduce the values of k ( i r ~ l) , i r = 1, n - r + 1 , at step p r in the values of 
k£ 2 \ = l,n-r + 2 

from step (p r _ j) . 

It follows: 



Jr-2) 






(r-1) 



fr 



7.(0 p-') fV) 7.1 r, _ (r-l) 7. W) 7.U-J 

K 1 T—> K n-r ’•••■> J„-r + 1 7 — <5/. 7 



Jr) 



7(0 
1 n-r + 1 



W) 



(0 



,('■-!) 7.(0 






/ r l = l,n - r - 1 , from which it follows that g ( '' 11 are linear functions with integer 
number coefficients. 

Then follows those (p r _ 2 ) in (p,._ e ) and so on, until we introduce the values 
obtained at step (p 2 ) in those from the step ( p t ) . 

It will follow: 




notation g. k x ,...,k n _ r , i =1 ,n , with all g. most obviously, linear functions with 

integer number coefficients (the notation was made for simplicity and aesthetical 
aspects). This is, thus, the general integer solution, of the initial system. 



The correctness of algorithm 1. 

The algorithm is finite because it has r steps on the forward way and r-1 steps 
on the reverse, (r < +oo) . Obviously, if one equation of one system does not have (integer 
number) solutions then the system does not have solutions either. 

Writing 5' for the initial system and S j the system resulted from step ( /r ) , 

\<j<r-2, it follows that passing from {pj) to (p J+l ) we pass from a system Sj to a 
system S j+l equivalent from the point of view of the integer number solution, i.e. 




ij = l,n-j + l. 



which is a particular integer solution of the system S / if and only if 



kl J) = h- 

7+1 7+1 



V+1 



= 1 ,n-j. 



is a particular integer solution of the system S J+l where 



L 0 _ r(j+\) 

' 7+1 ^ l j+\ 



1 ), i J+l = l,n-j. 



Hence, their general integer solutions are also equivalent (considering these 
substitutions). Such that, in the end, resolving the initial system S is equivalent with 
solving the equation (of the system consisting of one equation) 5) , with integer number 
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coefficients. It follows that the system S has integer number solution if and only if the 
systems S. have integer number solution, 1 < j < r — 1 . 



Example 1. By means of algorithm 1, let us calculate the integer number solution 
of the following system: 

f 5x - ly- 2z + 6w = 6 
(S) \ 

|-4.t + 6 v - 3z + 1 \w = 0 



Solution: We solve the first integer number equation. We obtain the general 
solution (see [4] or [5]): 



(A) 



x = t l + 2 1 2 

y = t i 

< 

z = -tj + 5i 2 + 3i 3 - 3 
w = t 3 



where t 2 , t, eZ. 

Substituting in the second, we obtain the system: 

(S { ) 5t 1 -23t 2 +2t 3 +9 = 0. 

Solving this integer equation we obtain its general integer solution: 



(Pi) 



h = K 

t 2 = 9A: t +23 k 2 + 7 



where k r k 2 eZ. 

The reverse way. Substituting ( p 2 ) in ( p ] ) we obtain: 



x = 3^ + 4^ 2 + 2 

y = K 

z = 31&j + 19k 2 +23 
w -9k x + 23k 2 + 7 

where k x , k 2 eZ, which is the general integer solution of the initial system (S') . Stop. 



Algorithm 2. 

Input 

A linear system (1) without all a {j = 0 . 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 

Method 

1. t = 1, h = 1, p = 1 

2. (A) Divide each equation by the largest co-divisor of the coefficients of the 
unknown variables. If you do not obtain an integer quotient for at least one equation, then 
the system does not have integer solutions. Stop. 
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(B) If there is an inequality in the system, then the system does not have integer 
solutions. Stop. 

(C) If repetition of more equations occurs, keep one and if an equation is an 
identity, remove it from the system. 

3. If there is (i 0 ,j 0 ) such that | a ioJo | = 1 then obtain the value of the variable 
x h from equation z 0 ; statement (T t ) . 

Substitute this statement (where possible) in the other equations of the system and 
in the statement (T t , ) , (H h ) and (P p ) for all /', h , and p . Consider t := t + 1 , remove 

equation (z 0 ) from the system. If there is no such a pair, go to step 5. 

4. Does the system (left) have at least one unknown variable? If it does, consider 
the new data and go on to step 2. If it does not, write the general integer solution of the 
system substituting k x , k 2 ,... for all variables from the right term of each expression 
which gives the value of the unknowns of the initial system. Stop. 

5. Calculate 

<7 = min \r\a.. =r mode/,.,. , 0 < Irl < Lx, 

j : : II y 1 y 2 II IJ2 

WW2 

and determine the indices as well as the r for which this minimum can be 

calculated. (If there are more variables, choose one arbitrarily.) 

Cl — T 

6. Write: x, = t h — x tj , statement ( H h ) . Substitute this statement (where 

Jl a 12 

S/2 

possible in all the equations of the system and in the statements ( T t ), ( H h ) and (P ) for 
all t, h , and p . 

7. (A) If a ^ 1 , consider x h := t h , h:=h + 1 , and go on to step 2. 

(B) If a = l, then obtain the value of x r from the equation (/) ; statement 
(P p ) • Substitute this statement (where possible in the other equations of the 
system and in the relations (T t ), ( H h ) and (P p _ x ) for all t, h , and p . 

Remove the equation (z) from the system. 

Consider h := h + 1, p := p + 1 , and go back to step 4. 

The correctness of algorithm 2. Let consider system (1). 

Lemma 1. We consider the algorithm at step 5. Also, let 

M= |r|, a ih = r mod a ih , 0 < \r\ < \a ih |, i, j v j 2 =1,2,3,... . 

Then M ^ 0 . 

Proof: 

Obviously, M is finite and M cz N* . Then, M has a minimum if and only if 
M ^ 0 . We suppose, conversely, that M = 0 . Then 
% = 0(moda (J . 2 ), V i,j\,j 2 . 

It follows as well that 

a i]2 = 0(moda l7| ), V 
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That is 



yiji | _ \ a ih |’ ^ A’ j 2 ' 

We consider an i 0 arbitrary but fixed. It is clear that 

: a i 0j * o. y/ 

(because the algorithm has passed through the sub-steps 2(B) and 2(C). But, because it 
has also passed through step 3, it follows that 

V/’ 

but as it previously passed through step 2(A), it would result that 

M = l > v/' • 

This contradiction shows that the assumption is false. 

Lemma 2. Let’s consider cl , = r(mod cl. ) . Substitute 

l oJ\ x v 2 



in system (A) obtaining system (B). Then 

o • ; 

Xj=Xj,j = \,n 

is the particular integer solution of (A) if and only if 

Xj = x%j * j 2 and t h = x° h - — — - 

a ioh 

is the particular integer solution of (B). 

Lemma 3. Let a x ^ and a, obtained at step 5. 

Then 0 < a 2 < a x 
Proof: 

It is sufficient to show that a x < |a, y |, V i, j because in order to get a 2 , step 6 is 
obligatory, when the coefficients if the new system are calculated, a, being equal to a 
coefficient form the new system (equality of modules), the coefficient on (i 0 j ] ) . 

Let a. . with the property \a, , \< a, . 

l oJo r r j | i q Jq | t 

Hence, a, > a, , = min \a. . . Let a, with a, , > a,, ; there is such an element 

1 | | \ l oJ\ l 0 Js I l oJs I I l Jm I ’ 

because a, , is the minimum of the coefficients in the module and not all W , , / = 1 ,n 

are equal (conversely, it would result that ) □ a jp V/el, r, the algorithm 

passed through sub-step 2(A) has simplified each equation by the maximal co-divisor of 
its coefficients; hence, it would follow that |a i ; | = 1, V/ = 1 ,n , which, again, cannot be 
real because the algorithm also passed through step 3). Out of the coefficients a,. . we 
choose one with the property a. igjs ^ Ma kj there is such an element (contrary, it would 

result (a j j,...,a ion )n | a ioJ | but the algorithm has also passed through step 2(A) and it 
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would mean that | a inj | = 1 which contradicts step 3 through which the algorithm has also 
passed). 



Considering q 0 = 



a / a . . 

‘0J So ‘0 Jm 



and 



r = a, . 

' oJs 0 



q 0 a 



’0 Jm 



we have 



Cl, 



l oJs 0 



s/- 0 (modfl. / ) and 0<|r 0 |< a inL < a kh < a x . We have, thus, obtained an r 0 with 



r 0 | < a x , which is in contradiction with the very definition of a x . Thus a x < a Xj , V/, / . 



Lemma 4. Algorithm 2 is finite. 

Proof: 

The functioning of the algorithm is meant to transform a linear system of m 
equations and n unknowns into one of m x x n x with m x < m , n x < n , thus, successively 
into a final linear equation with n—r + 1 unknowns (where r is the rank of the 
associated matrix). This equation is solved by means of the same algorithm (which works 
as [5]). The general integer solution of the system will depend on the n - 1 integer 
number independent parameters (see [6] - similar properties can be established also the 
general integer solution of the linear system). The reduction of equations occurs at steps 
2, 3 and sub-step 7(B). Step 2 and 3 are obvious and, hence, trivial; they can reduce the 
equation of the system (or even put an end to it) but only under particular conditions. The 
most important case finds its solution at step 7(B), which always reduces one equation of 
the system. As the number of equations is finite we come to solve a single integer number 
equation. We also have to show that the transfer from one system m i x n i to another 
m j+l x n i+l is made in a finite interval of time: by steps 5 and 6 permanent substitution of 
variables are made until we to a = 1 (we to a = 1 because, according to lemma 3, all 
a — s are positive integer numbers and form a strictly decreasing row). 

Theorem of correctness. 

Algorithm 2 correctly calculates the general integer solution of the linear system. 

Proof: 

Algorithm 2 is finite according to lemma 4. Steps 2 and 3 are obvious (see also 
[4], [5]). Their part is to simplify the calculations as much as possible. Step 4 tests the 
finality of the algorithm; the substitution with the parameters k x , /c,,... has 
systematization and aesthetic reasons. The variables t, h,p are counter variables (started 
at step 1) and they are meant to count the statement of the type T,H,P (numbering 
required by the substitutions at steps 3, 6 and sub-step 7(B); h also counts the new 
(auxiliary) variables introduced in the moment of decomposition of the system. The 
substitution from step 6 does not affect the general integer solution of the system (it 
follows from lemma 2). Lemma 1 shows that at step 5 there is always a , because 
0*McN\ 

The algorithm performs the transformation of a system m i x n j into another 
m i+x x n i+x , equivalent to it, preserving the general solution (taking into account, 
however, the substitutions) (see also lemma 2). 

Example 2. Calculate the integer solution of: 
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\2x-ly + 9z = 12 
- 5y + 8z + 10 w = 0 
0 z+ Ow = 0 
15x + 21z + 69w=3 

Solution: 

We apply algorithm 2 (we purposely selected an example to be passed through all 
the steps of this algorithm): 

1. t = 1, h = 1, p = 1 

2. (A) The fourth equation becomes 5x + 7z + 23w=l 

(B) - 

(C) Equation 3 is removed. 

3. No; go on to step5. 

5. a = 2 and i = 1, j\ = 2, j 2 = 3 , and r = 2 . 

6. z = +y, the statement (//,). Substituting it in the 

12x-2y+9q =12 

< 3y + 9q+10w =0 
5x +ly + 7t l +23w =1 

7.o^l consider z = t l ,h:=2, and go back to step 2. 

2.- 

3. No. Step 5. 

5. a = 1 and i = 2, y, = 4, j 2 = 2 , and r = 1 . 

6. y = t 2 -3w, the statement (7/ 2 ) . Substituting in the system: 

-12x + 2 1 2 + 9t y -6w = 12 

< 3 1 2 +8 1, + w = 0 
5x +7t 2 +7q +2w = 1 

Substituting it in statement (H^, we obtain: 
z = t x + 1 2 — 3w , statement (7/,) . 

7. w = -37,-8/) statement (7)). 

Substituting it in the system, we obtain: 

J-12x-20r 2 +57/j =12 
| 5x + 1 2 - 9t { = 1 

Substituting it in the other statements, we obtain: 
z = 10/, + 25/), (//,)" 
y = 10t, +24/p (//,)" 
h := 3, p\=2, and go back to step 4. 

4. Yes. 

2.- 

3. t 2 = 1 - 5x + 9fj , statement (T,) . 

Substituting it (where possible) we obtain: 
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{-1 12x + 237 t x = -8 (the new system); 

z = 10-50x + 115*p (//,)"' 
y= 10-50x + 114* p (H 2 ) n 
w = -3 + 15x + 35c , ( P 1 )' 

Consider t := 2 go on to step 4. 

4. Yes. Go back to step 2. (From now on algorithm 2 works similarly with that 
from [5], with the only difference that the substitution must also be made in the 
statements obtained up to this point). 

2 .- 

3. No. Go on to step 5. 

5. a = 13 (one three) and i = 1, j\ = 2, j 2 = 1 , and r = 13 . 

6. x = * 3 + 2t x , statement (77,) . 

After substituting we obtain: 

{- 1 12 C + = (the system) 

z = 10-50*3 +15f p 
y = 10-50C+ 14c, W"; 
w = -3 + 15C -5* p (/>)"; 

* 2 = 1- 5*3 - C ( 7 ) ) 

7. x := * 3 , /? := 4 and go on to step 2. 

2 .- 

3. No. go on to step 5. 

5. a = 5 and * = 1, /, = 1, j 2 = 2 and r = 5 

6. *, = * 4 + 9 * 3 , statement (H 4 ) . 

Substituting it, we obtain : 

5*3 + 13* 4 = -8 (the system). 



Z— 10 + 85*3+15*4, 


(tf.f; 


y = 10 + 76*3 +14* 4 , 


(// 2 ) /v 


x = 19*3+ 2* 4 , 




w = -3-30*3 ~ 3* 4 , 


(^i)'"; 


* 2 = 1-14*3 - C 


(7))"; 


:= 5 and go back to step 2. 



2 .- 

3. No. Step 5. 

5.o = 2 and i = 1, j\ = 2, j 2 = 1 and r = -2 . 

6. * 3 = *- - 3* 4 statement (H s ) . After substituting, we obtain: 
5* 5 -2* 4 = -8 (the system). 
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z = 10 + 85f 5 -240 1 4 , 


(vy 


y = 10 + 76t 5 -214r 4 , 


(H 2 ) v : 


x= 19t 5 -55t 4 , 


(H 3 ) iv : 


w = -3 - 30 1 5 +85t 4 , 


(py-. 


i 2 = — 1 — 1 4/, + 41t 4 , 


(7i)"'; 


t 4 = 9t 5 +26t 4 , 


(77 4 )'; 



7. t 3 := t 6 ,h := 6 and go back to step 2. 

2.- 

3. No. Step 5. 

5.o = l and i = 1, j { = 2, j 2 ,r = 1 . 

6. =t 6 + 2 1 5 statement (H 6 ) . After substituting, we obtain: 

t 5 - 2 1 6 = -8 (the system) 
z = 10-395t 5 -240t 6 , {H X ) VII \ 
y = 10-392t 5 -214t 6 , (T/,) 71 '; 
x = -91t 5 -55t 6 , ( H 3 ) ; 

w = -3 + 140t 5 + 85t 6 , (Pjf; 
f 2 = i + 68r 5 +41t 6 , (T,)" 7 ; 

t, = -43t 5 ~26t 6 , (H 4 )"; 

t 3 = -5t 5 -3t 6 , (H 5 ); 

7. t 5 =2t 6 — 8 statement (73) • Substituting it in the system we obtain: 0=0. 
Substituting it in the other statements, it follows: 

z = -1030t 6 + 3170 
y = -918f 6 + 2826 
x = -237 1 6 + 728 
w= 365t 6 -1123 
t 2 = 111 t 6 -543 
t 4 = 112t 6 +344 

> statements of no importance. 

t 3 = 13t 6 + 40 
t 4 = 5 1 6 - 16 

Consider h := l,p := 3 , and go back to step 4. t 6 eZ 

4. No. The general integer solution of the system is: 
x = -2377, +728 
y = -918kj +2826 
z = 1030k, +3170 
w = 365^ —1123 
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where k x is an integer number parameter. 

Stop. 

Algorithm 3. 

Input 

A linear system (1) 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 

Method 



1. Solve the system in □ ” . If it does not have solutions in □ " , it does not have 
solutions in Z" either. Stop. 

2. f=\,t = \,h = \,g=\ 

3. Write the value of each main variable x t under the form: 

(E f , ) : x, = + <h + 7 A i 



with all q tj , q t , r ip r n A,, in Z such that all |/v| <|A,|, A, ^ 0, |?*| <|A,| (where all x'j of 

the right term are integer number variables: either of the secondary variables of the 
system or other new variables introduced with the algorithm). For all i , we write 



4. 






r y x j 



njf.i+n 



= 0 where fed are auxiliary integer number 



variables. We remove all the equations [F f . J which are identities. 

5. Does at least one equation (F f . ) exists? If it does not, write the general 

integer solution of the system substituting k x ,k 2 ,... for all the variables from the right 
term of each expression representing the value of the initial unknowns of the system. 
Stop. 

6. (A) Divide each equation (f f by the maximal co-divisor of the coefficients 

of their unknowns. If the quotient is not an integer number for at least one i 0 the system 
does not have integer solutions. Stop. 

(B) Simplify -as in m - all the fractions from the statements (*/,<)■ 

7. Does r t h exist having the absolute value 1? If it does not, go on to step 8. If it 



does, find the value of x ; . from the equation \F f ^ J; write (D for this statement, and 
substitute it (where it is possible) in the statements (£/.,) ( r "‘) (ff.) (g«) for all 

i, t, h and g . Remove the equation (F f i(> Consider / := / + 1 , t := t + 1 , and go back 

to step3. 

8. Calculate 
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a- min r, r.. =r(modr ; ), 0< r < r.. 

; j ; II Vl l J2 II V2 

and determine the indices i m , j x , j 2 as well as the r for which this minimum can be 
obtained. (When there are more variables, choose only one). 

a i ~r , 

9. (A) Write x, = z h — x, , where z, is a new integer variable; statement 

a 

JmJl 

(B) Substitute the letter (where possible) in the statements 
( E f , , ) ( F f . , ) ( T , ) (H k _ j), (G g ) for all i, t, h and g . 

(C) Consider /?:=/? + 1 . 

10. (A) If a ^ 1 go back to step 4. 

(B) If a = 1 calculate the value of the variable x) from the equation 
relation (g', Substitute it (where possible) in the statements ( £ /,<) (A (<v,) 

for all i, t, h , and g . Remove the equation (f f . ). Consider g := g+ 1, / := / + 1 and go 
back to step 3. 



The correctness of algorithm 3 



Lemma 5. Let i be fixed. Then 



Z'rC +r 



V ■/=«! 



A ; (with all r., r t , A ; , « 2 being 



integers, « 1 < n 2 , A, ^ 0 and all x'. being integer variables) can have integer values if 
and only if r,, 2 , A, I ?;. 

Proof: 

The fraction from the lemma can have integer values if and only if there is a 



zeZ such that 



j 1 
\j='h j 



J ]r i jXj+r j I A, = z O £ r jjXj - A jZ + r t =0, 



j= n 1 



which is a linear equation. This equation has integer solution <=> r u ,..., r Uu , A, I r . 



Lemma 6. The algorithm is finite. It is true. The algorithm can stop at steps 1,5 or 
sub- steps 6(A). (It rarely stops at step 1). 

One equation after another are gradually eliminated at step 4 and especially 7 and 
10(B) (F f , ) - the number of equation is finite. 

If at steps 4 and 7 the elimination of equations may occur only in special cases 
elimination of equations at sub step 10 (B) is always true because, through steps 8 and 
9 we get to a = 1 (see [5]) or even lemma 4 (from the correctness of algorithm 2). 
Hence, the algorithm is finite. 

Theorem of Correctness. 

The algorithm 3 correctly calculates the general integer solution of the system (1). 
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Proof: 

The algorithm if finite according to lemma 6. It is obvious that the system does 
not have solution in □ " it does not have in Z" either, because Z" <z □ " (step 1). 

The variables /, t, h, g are counter variables and are meant to number the 
statements of the type E, F, T, H and G , respectively. They are used to distinguish 
between the statements and make the necessary substitutions (step 2). 

Step 3 is obvious. All coefficients of the unknowns being integers, each main 
variable x. will be written: 

( \ 



X; = 






C..X.+C,. 



I A, 



V i J 

which can assume the form and conditions required in this step. 

Step 4 is obtained from 3 by writing each fraction equal to an integer variable 
Y f i (this being x t eZ). 



Step 5 is very close to the end. If there is no fraction among all \E f i j it means 



that all main variables x ; already have values in Z , while the secondary variables of 
the system can be arbitrary in Z , or can be obtained from the statements T, H or G 
(but these have only integer expressions because of their definition and because only 
integer substitutions are made). The second assertion of this step is meant to 
systematize the parameters and renumber; it could be left out but aesthetic reasons 
dictate its presence. According to lemma 5 the step 6(A) is correct. (If a fraction 
depending on certain parameters (integer variables) cannot have values in Z , then the 
main variable which has in the value of its expression such a fraction cannot have 
values in Z either; hence, the system does not have integer system). This 6(A) also 
has a simplifying role. The correctness of step 7, trivial as it is, also results from [4] 
and the steps 8-10 from [5] or even from algorithm 2 (lemma 4). 

Ther initial system is equivalent to the “system” from step 3 (in fact, (£/,,) as 

well, can be considered a system) Therefore, the general integer solution is preserved 
(the changes of variables do not prejudice it (see [4], [5], and also lemma 2 from the 
correctness of algorithm 2)). From a system m i x n j we form an equivalent system 
m i+l x n i+l with m M < m i and n i+l < n i . This algorithm works similarly to algorithm 
2 . 



Example 3. Employing algorithm 3, find an integer solution of the following 
system: 



f 3xj + 4x, + 22x 4 - 8x 5 = 25 



< 6x x + 46x 4 - 12x 5 = 2 

4x 2 + 3x 3 - x 4 + 9x 5 = 26 

Solution 

1. Common resolving in □ 3 it follows: 
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*1 = 



x 3 = 



23x 4 - 6x s - 1 
-3 

x 4 + 2x 3 + 24 
4 

1 lx 5 + 2 
3 



2. f=l,t = l,h = l,g = l 



3. 



1 

H 

(N 




x, = lx 4 + 2x 3 + 

-3 


(+.) 


„ x, + 3x. 




x 2 = 6+ 4 5 

4 


(+0 


x s + 2 




X 3 - 4x 5 + ^ 


fr. 3 ) 



4. 



2* 4 +3y n -1 = 0 F l , 

■ * 4 + 2x 5 -4y 12 =0 F 1>2 

x 5 -3_v 13 +2=0 Fj 3 



5. Yes. 

6 . - 

7. Yes: |r 35 | = 1 . Then x 5 = 3y 13 -2 the statement (T t ) . Substituting it in the 
others, we obtain: 





_ , . 2x 4 - 1 

x, = 7x 4 + 6_y 13 4 + 


(+.) 






< 


*,= 6 + ^ +6 ^- 4 

4 


(+3) 




3v,, — 2 + 2 

x 3 = -12 v 13 + 8 + -^ 


<+ 0 


Remove the equation (FJ 3 J. 




Consider / := 2 


, 2 ; go back to step 3. 






_ , . 2x 4 - 1 

x, = 7x 4 + 6_y 13 4 + 


(+.) 


3 




(+ 




x 3 = - 1 l+ l3 + 8 


(E 2 , ,) 
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4. 



2*4 +3.V 21 -1 = 0 F 2 , 

< 

*4 + 2y,3 ~ 4y 22 = 0 F 2 2 

5. Yes. 

6 . - 

7. Yes |r 24 | = 1 . We obtain x 4 = -2_y l3 + 4y 22 statement (T 2 ) . Substituting it in 
the others we obtain: 



v — ■ 


-28y 22 + 20y 13 + - 


-4*3+8^22-1 


F 


*1 


-3 


^ 2 , 1 


1 = 


y 22 + * 3 + 5 




F 

^ 2 , 2 


*3 = 


-11*3+8 




F 

^ 2 , 3 



Remove the equation (F 2 2 ) 

Consider / := 3, t := 3 and go back to step 3. 

3. 

*, = -22y„ + 30y„ + 2 - V|1+ _ 2 3 > ' i2 - 1 E, , 

*2 = *3 + -V 2 2 + ^ E 2 * 2 

x 3 = - 1 ly l3 +8 £ 33 

4. 2y 13 + 2y 22 + 3y 31 -l=0 F xl 

5. Yes. 

6 . - 

7. No. 

8. a = l and i m = 1, /, =31, j 2 = 22 , and r = 1 . 

9. (A) _y 9 , =Zj — y 3l (statement (// )). 

(B) Substituting it in the others we obtain: 



Xj = -22y 13 - 30 z 3 + 30y 31 


/j | 2y l3 +2 Zl -2y 3l -l 
-3 


E 

^3, 1 


*2 = *3+ Zl “* 


+ 5 


E x 2 ' 


*3 “ 1 !* 3 


+ 8 


^ 3 . 3 ) 


2*3 + 2zi + y 31 -l = 0 


^3,1 ’ 




*4=-2* 3 +4Zj-4y 13 


t 2 ' 





(C) Consider h := 2 

10. (B) y u = 1 -2y u -2 z l , statement (G 3 ). 
Substituting it in the others we obtain: 
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x, = -40y, 3 -92zj+27 


F 

^3,1 


x 2 = 3y, 3 + 3z, + 4 


F 

^3,2 


u* 

II 

1 

> r 

u> 

+ 

OO 


F 

^3,3 


x 4 = 6y 13 + 12zj -4 


T 2 


^22 =2.Vi3+3Zi -1 





Remove equation (F 3 , ). 

Consider g := 2, / := 4 and go back to step 3. 
3. 



x, = -40y, 3 - 92 z, + 27 


^4,1 


x 2 = 3y 13 + 3z, + 4 


^4,2 


& 

ii 

i 

U) 

+ 

OO 


^4,3 



4. - 

5. No. The general solution of the initial system is: 



x, = -40k, - 92k 2 + 27, 


from 


£ 4 ,i 


Xt = 3k, + 3 k 2 +4, 


from 


^4,2 


li 

i 

+ 

OO 


from 


F 

^4,3 


x 4 = 6k, +12k 2 -4, 


from 


r T~' ft 

L 2 


x 5 = 3k, -2, 


from 


0 



where k,,k 2 eZ. 

Algorithm 4 
Input 

A linear system (1) with not all a tj = 0 . 

Output 

We decide on the possibility of integer solution of this system. If it is possible, we 
obtain its general integer solution. 

Method 

1. h — 1, v = 1 . 

2. (A) Divide every equation i by the largest co-divisor of the coefficients of the 
unknowns. If the quotient is not an integer for at least one i 0 then the system 
does not have integer solutions. Stop. 

(B) If there is an inequality in the system, then it does not have integer 
solutions 

(C) In case of repetition, retain only one equation of that kind. 
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(D) Remove all the equations which are identities. 

3. Calculate a = min Lx. , a u ^0 and determine the indices i 0 ,j 0 for which this 

i,j i 1 

minimum can be obtained. (If there are more variables, choose one, at 
random.) 

4. If a ^ 1 go on to step 6. 

If a - 1 , then: 

(A) Calculate the value of the variable x h from the equation i 0 note this 
statement C> 

(B) Substitute this statement (where possible) in all the equations of the 
system as well as in the statements (V vl ), (// ;i ), for all v and h . 

(C) Remove the equation i 0 from the system. 



(D) Consider v := v + 1 . 

5. Does at least one equation exist in the system? 

(A) If it does not, write the general integer solution of the system substituting 
k x , k 2 ,... for all the variables from the right term of each expression 
representing the value of the initial unknowns of the system. 

(B) If it does, considering the new data, go back to step 2. 

6. Write all <?,. /' ^ /' and b ; under the form : 

l oJ J J u *0 



a, = a , , tr , + r . , with r. < a 

l oJ l oJo ll 0 j *0 j I in / I I 

b, . = a , q, + r , with I r t I <| a, 

*n./n ± ‘o 'n ‘n 



*0 J 



l oJo 



7. Write x h = - £ 7, 0 A ; + q io + h > statement (/7, ). 

i*i o 

Substitute (where possible) this statement in all the equations of the system as 
well as in the statement (// /( ), for all v and h . 

8. Consider 

x jo := t h , h := h + 1, 



hJ 



* j 0 > 



a t := ±a. . 

*070 *070 



K '■= +n, ^ 



and go back to step 2 



The correctness of Algorithm 4 

This algorithm extends the algorithm from [4] (integer solutions of equations to 
integer solutions of linear systems). The algorithm was thoroughly proved in our previous 
article; the present one introduces a new cycle - having as cycling variable the number of 
equations of system - the rest remaining unchanged, hence, the correctness of algorithm 
4 is obvious. 
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Discussion 

1. The counter variables h and v count the statements H and V , respectively, 
differentiating them (to enable the substitutions); 

2. Step 2 ((A)+(B) + (C)) is trivial and is meant to simplify the calculations (as 
algorithm 2); 

3. Sub-step 5 (A) has aesthetic function (as all the algorithms described). 
Everything else has been proved in the previous chapters (see [4], [5], and 
algorithm 2). 

Example 4. Let us use algorithm 4 to calculate the integer solution of the 
following linear system: 

f3xj -7x 3 +6x 4 =— 2 

|4xj + 3i 2 + 6x 4 -5x 5 =19 

Solution 

1. h = 1, v = 1 

2 . - 

3. a = 3 and i = 1, j = 1 

4. 3 ^ 1 . Go on to step 6. 

6. Then, 

-7 = 3- (-3) + 2 
6 = 3 -2 + 0 
-2 = 3 • 0 - 2 

7. Xj =3 x 3 - 2x 4 + t l statement («,)• Substituting it in the second equation we 
obtain: 

4fj + 3 x 2 + 12x 3 - x 4 - 5x 5 =19 

8. Xj := t u h := 2, a l2 := 0, a n := +2, a 14 := 0, a n := +3, b :=—2 . 

Go back to step 2. 

2. The equivalent system was written: 

J3q +3 x 3 =— 2 

[4q +3 x 2 +12x 3 — x 4 -5x 5 =19 

3. a = 1, i = 2, j = 4 

4. 1=1 

(A) Then: x 4 = 4/, +3x 2 +12x 3 -5x 5 -19 statement (Ej). 

(B) Substituting it in (TEj), we obtain: 

Xj = — 7q - 6 x 2 - 2 lx 3 + 10x 5 + 38, (».) 

(C) Remove the second equation of the system. 

(D) Consider: v := 2 . 

5. Yes. Go back to step 2. 

2. The equation +3 q + 2x 3 = -2 is left. 

3. a = 2 and i = 1, / = 3 
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4. 2 + 2 , go to step 6. 

6. 

+3 = +2-2-1 
-2 = +2(-l) + 0 

7. +, = -2?j + t 2 — 1 statement (// 2 ). 

Substituting it in , (k,), we obtain: 

Xj = 35f, - 6 x 2 - 2 1? 2 + 10x 5 +59 H x " 
x 4 = -20f, +3x 2 + 12f 2 -5x 5 -31 V x ' 

8. x 3 := t 2 , h := 3, a n := — 1, a 13 := +2, 6^=0 , (the others being all = 0). Go back 
to step 2. 

2. The equation -5/\ +2t 2 =0 was obtained. 

3. o = l, and i = 1, / = 1 

4. 1=1 

(A) Then q =2^ statement (l^). 

(B) After substitution, we obtain: 

Xj = 49t 2 -6x 2 +10x 5 +59 H x 

x 4 = -28t 2 +3x 2 -5x 5 -31 V 7 " ; 
x 3 — — 3t 2 H 2 , 

(C) Remove the first equation from the system. 

(D) v := 3 

5. No. The general integer solution of the initial system is: 

X[ = 49 k x - 6k 2 + 10£ 3 + 59 
x 2 = k 2 

x 3 = -3&j - 1 

x 4 = -28£j + 3k 2 - 5k 3 - 3 1 
x 5 = £ 3 

where (/q, k 2 ,k 2 )eZ i . 

Stop. 

Algorithm 5 
Input 

A linear system (1) 

Output 

We decide on the possibility of an integer solution of this system. If it is 
possible, we obtain its general integer solution. 

Method 

1. We solve the common system in □ " , then it does not have solutions in □ " , 
then it does not have solutions in Z” either. Stop. 

2. f = l,v=l,h = l 

3. Write the value of each main variable x, under the form: 
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f 



E f.i '/,+ Zv v '/+' 



A 



J 



/A ; , 



7 V 1 

with all q tj , q t , r.., r t , A ( . from Z such that all |r.|<|A,.L I/; I < |A,| , A, ^ 



4 . 



5. 



6 . 



7. 

8 . 



9. 



(where all v ; - S of the right term are integer variables: either from the 

secondary variables of the system or the new variables introduced with the 
algorithm). For all i , we write r = A, 



(a)=Z 



r s x i 



'■•Jr 



y f . 



+ /;. = 0 where are auxiliary integer variables. 



Remove all the equations [F f J which are identities. 

Does at least one equation (iy , ) exist? If it does not, write the general integer 

solution of the system substituting k v k 2 ,... for all the variables of the right 
number of each expression representing the value of the initial unknowns of 
the system. Stop. 

(A) Divide each equation (F f . ) by the largest co-divisor of the coefficients of 



their unknowns. If the quotient is an integer for at least one i 0 then the system 
does not have integer solutions. Stop. 

(B) Simplify - as previously ((A)) all the functions in the relations (£/..)■ 



Calculate <7 = min Ir.l, r. ^0 , and determine the indices i 0 ,j 0 for which this 

i,j I 1 

minimum is obtained. 

If a ^ 1 , go on to step 9. 

If <7 = 1, then: 

(A) Calculate the value of the variable x h from the equation (A) write 
(V v ) for this statement. 

(B) Substitute this statement (where possible) in the statement (k f t 
(A). (h„), for all i, v , and h . 

(C) Remove the equation (E f . 

(D) Consider v := v + 1, / := /+ 1 and go back to step 3. 

Write all r . , 7 ^ /' and r. under the form: 

r, , = A, • q, + r , with \r . < I A, I ; 

l oJ l 0 1 l oJ l oJ | z 0-/ I I 1 ' 

r = A, • q, + r. , with r' < I A, I . 

l oJ l o 1 l o l o Mo I in 



10. (A) Write x' Jo = + % + h statement H h . 

j*h 

(B) Substitute this statement (where possible) in all the statements 

(a) ( k )' C.-,)' 
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(C) Consider h := h + 1 and go back to step 4. 

The correctness of the algorithm is obvious. It consists of the first part of 
algorithm 3 and the end part of algorithm 4. Then, steps 1-6 and their correctness were 
discussed in the case of algorithm 3. The situation is similar with steps 7-10. (After 
calculating the real solution in order to calculate the integer solution, we resorted to the 
procedure from 5 and algorithm 5 was obtained). This means that all these insertions 
were proven previously. 



Example 5 

Using algorithm 5, let us obtain the general integer solution of the system: 
f 3x, + 6x 3 + 2x 4 = 0 

i 4x 2 -2x 3 - lx 5 = -1 



Solution 

1. Solving in □ 5 we obtain: 



< 






x. 



-6x 3 - 2x 4 
_ 

-2x 3 + 7x 5 -1 

4 



2 . 

3. 



/'= 1 , v= 1 , h = \ 




: Xj = 2 x 3 + 



-2x 4 

3 




x 5 + 



2x 3 + 3x 5 - 1 
4 



4. F u :-2x 4 -3y M =0 

(^,2 ): 2x 3 + 3x 5 - 4 y 12 -1 = 0 

5. Yes 

6 . - 

7. i = 2 and i 0 = 2, j 0 = 3 

8 . 2*1 

9. 3 = 24 + 1 
-4 = 2 -(-2) 



-1 = 2 - 0-1 

10. x 3 =-x 5 + 2+, 2 +t t statement (//, ). After substitution: 



E h i 

E \'2 



' : x, = 2x 5 - 4 v 12 - 2t t + - 
' : x 2 = x 5 + 



-2x 4 

3 

x 5 +4y 12 + 2r, -1 

4 



F 12 ' : x 5 + 2r, - 1 = 0 
Consider h := 2 and go back to step 4. 
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4. F u ' : —2 x 4 —3y n =0 

F l2 ' \2t\+x 5 -\ = 0 

5. Yes. 

6 . - 

7. a = 1 and i 0 = 2, j 0 = 5 

(A) v 5 = -2 t l + 1 statement (V t ) 

(B) Substituting it, we obtain: 

E u x l =-6» l +2-4v l2 +^i- 

^ 1,2 ’ -^2 = ~ +1 + Ji 2 

7/j x 3 = 3fj +1-1 + 2 y 12 

(C) Remove the equation 2 ). 

(D) Consider v = 2, / = 2 and go back to step 3. 

3. : x,=-6t,-4y, 1 +2+^j± 

(^ 2,2 -^2 — — + kl 2 1 

4. F 2A : -2x 4 -3_y 12 =0 

5. Yes. 

6 . - 

7. a = 2 and z 0 = 1, j 0 = 4 

8 . 2*1 

9. - 3 = -2 • (1) - 1 

10. (A) x 4 = —y 21 +t 2 statement (//,) 

(B) After substitution, we obtain: 

E 2 , i ' : A = -6h - 4 v 12 + 2 + 2y ^ 

F 2 .i -y 2 i - 2 c =o 

Consider h := 3, and go back to step 4. 

4 - f 2,i -y 2 i-2c=o 

5. Yes 

6 . - 

7. a = l and z 0 = 1, j 0 = 21 (two, one). 

(A) y 2l =-2 1 2 statement (k,). 

(B) After substitution, we obtain: 

(C) Remove the equation F 2 , . 

(D) Consider v = 3, / =3 and go back to step 3. 
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3. (E xl ): ^ = -6 fj - 4 j 12 -2t 2 +2 

(^3,2 )• "U = — 2^ +JF 12 +1 

4. - 

5. No. The general integer solution of the system is: 



<N 

-U; 

1 

1 

II 


— 2 k 2 + 2 , 


from 


F 

^ 3,1 ’ 


x n — — 2k| + k 2 


+ 1, 


from 


F 

^ 3,2 ’ 


x, = 3 k x + 2 k 2 


- 1 , 


from 




x 4 = 


3 k 3 , 


from 


h 2 


x 5 = -2 k x 


+ 1 , 


from 


O 



where (k x , k 2 ,k^)e Z . 

Stop. 

Note 1. Algorithm 3, 4, and 5 can be applied in the calculation of the integer 
solution of a linear equation. 

Note 2. The algorithms, because of their form, are easily adapted to a computer 
program. 

Note 3. It is up to the reader to decide on which algorithm to use. Good luck! 
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